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Abstract 

Intermediate Mass Black Holes (IMBHs) are candidates to seed the Supermassive 
Black Holes (SMBHs), and some could still wander in the Galaxy. In the context of 
annihilating dark matter (DM), they are expected to drive huge annihilation rates, 
and could therefore significantly enhance the primary cosmic rays (CRs) expected 
from annihilation of the DM of the Galactic halo. In this proceeding (the original 
paper is Brun et al. 2007), we briefly explain the method to derive estimates of 
such exotic contributions to the p and e + CR spectra, and the associated statistical 
uncertainties connected to the properties of IMBHs. We find boost factors of order 
10 to the exotic fluxes, but associated with very large statistical uncertainties. 
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As in other topics in fundamental physics, annihilating dark matter can be 
motivated by some coincidence argument, because exotic microscopic physics 
is expected to surge at the electroweak scale: This energy scale gives ideal 
masses as well as good interaction strengths to naturally get the correct abun- 
dance of WIMPs today, provided that they are thermally produced in the 
early U niverse, without m atter-antimatter asym metry (see the ver y nice lec- 



ture by iMurayamal 120071 ) . In the original paper (IBrun et all 120071 ). we have 



studied the possibility that IMBHs could be the sources of a copious amount 
of DM annihilatio n products, more pre cisely p's and e + 's. Such sources have 



been proposed by iBertone et al.l (120051 ) in the fra me of ind i rect d etection of 
DM with 7-rays, and f urther extended to ne utrinos (iBertond . 120061 ) . Using the 
method developed by lLavalle et al.l (j2007al ). we have computed how IMBHs 
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may boost the exotic production of e +/ s and p's, compared to considering 
only an nihilation of the DM of th e Galactic smooth halo (study complemen- 
tary to iBringmann fc Salatill2007l ). In this proceeding, we shortly present the 
method to estimate the boost factors to the fluxes for both species and the 
corresponding variances associated with the statistical properties of IMBHs. 
The exhaustive study with some predictions for benchmark WIMP models 
(super symmetry or extra-di mensions) are a vailab le in the original paper. We 
concentrate on scenario B of iBertone et al.l (120051 ) for the formation of IMBHs 
during the matter era. In this scenario, IMBHs form out of cold gas in early- 
collapsed DM halos and are typified by a mass scale of order 1O 6 M . During 
their formation, if adiabatic, IMBHs modify the DM density profiles of their 
host halos, so that spikes appear at the centres. Starting with an r~ 7 spherical 
density profile, the final logarithmic slope due to the adiabatic compression 
of dark matter is 7f} n = (9 — 27)/ (4 — 7). For instance, taking a profile ini- 
tially with 7 = 1 leads to 7g n = 7/3, which stands inside the spike radius 
r sp . Whatever the DM profile surrounding an IMBH, it is useful to define 
a quantity proportional to the corresponding WIMP annihilation rate, the 

effective annihilation volume, as follows: £ = / BH d 3 x (y 11 ) , where pbh is 
the DM density profile around the black hole, and, p@, the local dark mat- 
ter density, allows for a normalisation to the local annihilation rate. Such an 
effective volume is the one the DM stuck to the IMBH would have if it were 
diluted down to the local density p & . Note that £ is sensitive to the inner 
cut-off radius r c below which the spike density saturates, which is usually set 
by equating the gravitational infall rate with the WIMP annihilation rate. 
This cut-off is associated with the maximum value of the accreted DM den- 
sity p max = p(r c ) ~ Tn x /((a ann v)r), where r ~ 10 Gyr is the typical IMBH 
formation timescale. If computed inside r sp , £ oc (m x /(cr ann i;)) 5 / 7 , and only 
slightly depends on the WIM P model. Cons i derin g; IMBHs forming in early- 



collapsing halos (r 1 profiles), IBertone et al.l (120051 ) performed simulations to 



derive the evolution of their DM spikes in the Galaxy. From their results, they 
extracted the statistical properties of the surviving spikes, i.e. their spatial 
distribution in the Galaxy, and the surrounding DM profile properties. The 
total number of IMBHs is small, ~ 100 within a galactocentric radius of 500 
kpc. Nonetheless, the probability to find an IMBH increases when going to- 
ward the Galactic Centre, and the Earth is located at an interesting place from 
an observational point of view. The whole phase space of IMBH properties is 
thus portrayed by the spatial distribution dVyjdV of objects in the Galaxy 
and the distribution of their effective annihilation volumes £, dV^/d^, both ex- 
tracted from these simulations, that we used. For generic WIMPs, the typical 
mean value (£) ~ 10 6 kpc 3 , which is huge. The CR propagation modelling is 
a key ingredient for this kind of studies. Here, we adopt a slab diffusion zone, 
featured by its radial extension that we fix to -R s i a b — 30 kpc, and by its half- 
thickness L, 4 kpc here. CRs are confined within the slab, which translates 
to Dirichlet boundary conditions for the diffusion equation (the CR number 



2 



density vanishes on the borders). For the transport processes, we take a spa- 
tial independent diffusion coefficient K(E) = f3K lZ s (where 71 = pc/Ze is 
the rigidity) and a constant wind Konv directed outwards along the vertical 
axis. Such a configuration is rem iniscent from the medium set of parameters 



provided by lMaurin et al.ll200ll . that we choose here: Kq = 0.0112 kpc /Myr 



S = 0.7 and V conv = 12 km/s. One can easily write and solve the diffusion 
equations for both e +/ s and p's for this kind of geometry. For e + 's, the main 
processes that come into play are the energy losses (mainly inverse Compton 
diffusion off CMB or IR photons, and synchrotron radiation), and the diffusion 
on the magnetic turbulences. Disregarding the convection process, which is 
much less efficient than energy losses, one can express the typical propagation 

length for e+'s as: A d = \ 2K r E y [zf^ JJ where e = E/{E = 1 GeV}. 

Assuming an infinite 3D spherical diffusion zone (correct while Ad < L), the 
e + propagator is proportional to a Gaussian function of the source distance 
with a = Ad (< few kpc): Sources located farther than Ad will almost not 
contribute to the flux at the Earth. A^ being a decreasing function of the de- 
tected energy, the effective volume in which e +/ s propagate increases as they 
loose energy. The diffusion equation for p's must include spallation processes 
and wind convection that occur in the thin Galactic disc, so we can not use a 
simple spherical symmetry to derive a global expression. Moreover, the energy 
losses are negligible for this species, which modifies significantly the picture 
that we had for e + 's. Nevertheless, it is useful to write, as for e + 's, the typi- 
cal propagation length: Ad = W^-, where convection is assumed to dominate 

vconv 

over spallation in average, which is correct unless at sub-GeV energies. This 
is quite different from e + 's because this length is an increasing function of 
energy. The picture is therefore reversed, and the propagation volume is much 
larger at higher energy for p's (Ad reaches the size of the diffusion slab at 
energies of order 10-100 GeV). We now define a convenient Green function 
for any CR species, that encodes the injected spectrum induced by DM an- 
nihilation: G(E,x Q <- x) = J^dEs G(E,x Q «- E s ,x) x dN ™f s \ where 
dNcR/dEs is the injected spectrum at source (E = Eg for p's). The Galactic 
host halo is described by a smooth DM distribution p sm , so that the corre- 
sponding primary CR flux reads: (f) sm (E) = J halo d 3 x G(E,x Q *- x) (^) 
where v is the CR velocity and S = S(a a _ nn v)p Q 2 /(2m x 2 ) encodes the main 
WIMP propertied]. IMBHs can be considered as point-like sources, and the 
flux due to the i th object is merely: (j)^i(E) = x S x £ 4 x Q(E,x Q <— Xj). 
Then, we have to sum over the whole population. We can derive a statis- 
tical prediction by using the IMBH phase space information, either with 
MC simulations or with a semi-analytic method. The overall IMBH contribu- 



tion is: bh , tot (£) = ^E^£f=i<M£) - i-JdNN^l x J 



M t-^i=\ *sj=l yi;3\ ) 47r J w J-* dN A J U SS d £ 



X 



Jd 3 xG(E,x G <- x)^m = (A)(0 bh ) = jL S (N)(0(G} where the 1 st equality 



5=1/2 for Majorana particles, 1 otherwise. 
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is a mere counting (performed e.g. with a Monte Carlo — MC), and the 2 nd 
one takes directly the — normalised — pdfs into account (this limit is reached 
for an infinite number M of MC realisations). Those pdfs characterise the total 
number N of IMBHs wandering in the Galactic halo, £ and Q, the latter being 
spatially weighted with dVy jdV . The last equality gives the same quantities in 
terms of statistical mean values. This assumes no correlations between the con- 
sidered variables, which has been carefully checked from the simulation results. 
The calculation of the variance of CR fluxes originating from the whole IMBH 

2 / 2 2 22\ ^ 

sample a bh)tot is straightforward: = ^ + + + > 

where a x is the variance corresponding to any variable x. The boost factors 
for e + 's and p's are the ratios of the fluxes originating from a halo popu- 
lated with IMBHs to those calculated for the host smooth halo alone. As 
the Galaxy mass fraction carried by IMBHs is negligible, the effective boost 
is B(E) ~ 1 + 0bh,tot/0sm, which depends on energy. Estimates of the mean 
fluxes for e + 's and p's and associated statistical variances have been performed 
by using (i) a MC method (ii) a se mi-analytic calculation based on the method 
proposed by lLavalle et al.l (j2007al ). Both methods agree, and the average boost 
factors are found to be ~ 10 4 for both species, but with a very large rela- 
tive variance, going up to ~ 100% as the CR propagation scale decreases ; 
those results are very different from what is found in the co ntext of stan- 



dard D M dumpiness, for which the enhancement is negligible (lLavalle et al. 



2007bl ). The large statistical error found is connected with the small number 



of IMBHs predicted inside the diffusion volume set by the CR propagation 
scale, and more precisely with the distance of the closest object. Should it 
be close, which is statistically unlikely, it would strongly dominate the exotic 
contribution of antimatter, otherwise it would remain almost unobservable. 
The existing measurements of the e + and p fluxes already provide strong 
constraints that we will translate in limits on the gamma-ray and neutrino 
production from IMBHs in a forthcoming paper. 

Acknowledgements: The author is grateful to his collaborators, G. Bertone, 
P. Brun, P. Salati and R. Taillet for enlightening discussions during this work. 



References 

Bertone, G., 2006, Phys. Rev. D 73 (10), 103519— h 
Bertone, G., Zentner, A. R., Silk, J., 2005, Phys. Rev. D 72 (10), 103517. 
Bringmann, T., Salati, P., 2007, Phys. Rev. D 75 (8), 083006. 
Brun, P., Bertone, G., Lavalle, J., Salati, P., Taillet, R., 2007, Phys. Rev. D 
76 (8), 083506-+. 

Lavalle, J., Pochon, J., Salati, P., Taillet, R., 2007, A&A 462, 827-840. 
Lavalle, J., Yuan, Q., Maurin, D., Bi, X.-J., 2007, ArXiv e-prints 0709.3634. 
Maurin, D., Donato, F., Taillet, R., Salati, P., 2001, ApJ 555, 585-596. 
Murayama, H., 2007, ArXiv e-prints 0704.2276. 



4 



